BCS ansatz for superconductivity in the canonical ensemble 
and the Pauli exclusion principle 



G. Zhu\ M. Combescot^'^ and O. Betbeder-Matibet^ 
(1) Department of Physics, University of Illinois at Urbana- Champaign, 1110 W Green St, Urbana, IL, 61801 and 
(2) Institut des NanoSciences de Paris, Universite Pierre et Marie Curie, 
CNRS, Tour 22, 4 place Jussieu, 75005 Paris 
(Dated: December 2, 2011) 

The usual formulation of the BCS ansatz for superconductivity in the grand canonical ensemble 
makes the handling of the Pauli exclusion principle between paired electrons straightforward. It 
however tends to mask that many-body effects between Cooper pairs interacting through the reduced 
BCS potential are entirely controlled by this exclusion. To show it up, one has to work in the 
canonical ensemble. Pauli blocking between a fixed number of composite bosons is however known 
to be difficult to handle. To do it, we hero develop a commutator formalism for Cooper pairs, 
along the line we used for excitons. We then rederive, within the A'^-pair subspace, a few results of 
BCS superconductivity commonly derived in the grand canonical ensemble, to evidence their Pauli 
blocking origin. We end by discussing what should be called "Cooper pair wave function" . 

PACS numbers: 



I. INTRODUCTION 

A major breakthrough in the understanding of 
superconductivity is definitely due to the wave 
function ansatz proposed by Bardeen, Cooper, 
SchriefFer [1]. The great advantage of this ansatz 
is that it leads to results in agreement with ex- 
periments through calculations quite easy to per- 
form analytically. Its standard form in the grand 
canonical ensemble as a sum of states with differ- 
ent particle numbers, makes the Pauli exclusion 
principle between paired up and down spin elec- 
trons straightforward to handle. Calculations in 
the grand canonical ensemble however mask the 
key role played by the Pauli exclusion principle in 
superconductivity. Indeed, due to the very pecu- 
liar form of the reduced BCS potential — an up 
spin electron (k) can interact with a down spin 
electron (— k) only — Pauli blocking is the only 
way two correlated electron pairs can feci each 
other (see Fig. 1). This in particular explains 
why Cooper pairs can strongly overlap witho ut 
dissociating, by contrast with excitons which dis- 
sociate into an electron-hole plasma [2] through 
a Mott transition when overlap starts. To bet- 
ter understand the key role played by the Pauli 
exclusion principle in the physics of BCS super- 
conductors, it is necessary to stay in the canon- 
ical ensemble, with a pair number and a num- 
ber of states available for pairing fixed. However, 
to exactly handle Pauli blocking between a fixed 
number of paired fermions like excitons or Cooper 
pairs, is far from easy. There were, in the past, 
several discussions on the exact particle-conserving 
solution within the reduced BCS potential using 
Richardson-Gaudin procedure [3-9] and its differ- 
ence with the BCS ansatz [10-14]. Most discus- 
sions on the ground state however focus on recover- 
ing the correct energy or some physical quantities. 




FIG. 1: Two Cooper pairs made of the repeated inter- 
actions of an up spin electron k and a down spin elec- 
tron — k cannot interact by the reduced BCS potential 
because this would impose k'l = k2. The two up spin 
electrons would then have the same momentum which 
is impossible due to the Pauli exclusion principle. 



like the integrals of motion ([8]), instead of the 
wave function itself. Yet, the BCS ansatz wave 
function is strongly linked to the picture people 
commonly have of superconductivity. 

Diiring the last decade, we developed a many- 
body formalism appropriate to composite bosons 
made of fermion pairs[7, 15-17]. Up to now, we 
have extensively studied excitons and developed 
a formalism adapted to their many-body physics. 
Because the long-range Coulomb interaction be- 
tween electrons and holes leads to a Mott dissocia- 
tion of the exciton gas into an electron-hole plasma 
[2] when the density increases, the relevant exciton 
regime is the dilute regime. By contrast, the rele- 
vant regime in BCS superconductivity is the dense 
regime with Cooper pair wave functions strongly 
overlapping. As a result, the many-body physics 
of composite bosons like Cooper pairs is expected 
to have some similarities with the one of excitons, 
but with a few major differences. 

In section II, we develop a commutator formal- 
ism for paired electrons capable to handle the Pauli 
exclusion principle within a N-pah condensate, in 
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an exact way. 

In section III, we come back to the BCS ansatz 
in its grand canonical form and we briefly rcdcrive 
through the standard grand canonical procedure, 
some textbook results [18, 19] on quantities we 
are going to consider in the canonical ensemble, 
namely, the k-state population, the pair number 
mean value, its fluctuations, and the two-pair cor- 
relation function. 

In section IV, we turn to the canonical ensem- 
ble, with a fixed number of pairs. Through a di- 
rect study of the probability distribution for N- 
pair states in the BCS ansatz, we prove that this 
ansatz indeed corresponds to a distribution very 
much peaked on a particular value N* of the pair 
number, as a result of the "moth-eaten effect" 
induced on Cooper pairs by the Pauli exclusion 
principle [17]. The standard derivation of this re- 
sult, through the fluctuations of the pair number 
around its mean value[18], completely hides the 
microscopic origin of this maximum. 

In section V, we calculate the fraction of the 
k electron state occupied in a A''-pair state. We 
do show that it is identical to the one calculated 
within the grand canonical version of the ansatz 
for N = N* . This is also true for the pair operator 
mean value associated to what is often called "pair 
wave function" . 

In section VI, we come back to what should be 
called "Cooper pair wave function". 

In section VII, we conclude. 



II. COMPOSITE BOSON FORMALISM 
FOR CONDENSED PAIRS 

The goal of this section is to develop a formalism 
capable to, in an exact way, handle Pauli blocking 
within a A^-fermion-pair condensate. 

For that, we follow our previous works [7, 17] 
and introduce the generalized creation operator for 
correlated pairs. 



with n = (0, 1,2---). The operator = alj)^_i^ 
creates a pair of free fermions with zero total mo- 
mentum. In the case of BCS superconductivity, 
these fermions are up and down spin electrons. 
We first note that free fermion pair creation op- 



erators commute, [/3^, , pI.] = 0, while 
[/3k',/3l] = <5k',k--Dk'k , 



(2) 



with Dk'k = '^k',k(akak + b^_]J)-k), so that 
£>k'k|0) = 0. We also have 

[(iW, 4] = ^P,k4 = [blpb.p, Pi] ■ (3) 



It is then easy to show that 



(4) 



where the "deviation from boson operator" of these 
generalized correlated pairs is given by Dm = 
Ek kkr+Halok + bl^^b-k). The scalar t„, de- 
fined as 



m+l 



(5) 



is the m+l moment of the k state distribution 
in the correlated pair at hand. To possibly re- 
late this moment to the correlated pair wave; func- 
tion, we are led to normalize the ifk distribution 
as To = J2k Iv'kl ^ then have, for |0) being 

the vacuum state 



(0|Bo4lO)=To = l 



(6) 



In order to easily handle the Pauli exclusion 
principle within a Bq^|0) condensate, we also need 



[Dm,Bi]=2B^ 



m+n+1 • 



(7) 



Using it and 



n R' 



B, 



fJV-l 




+bI Dm,Bl 



we get by iteration 



Dniy Bq 



(8) 



(9) 



In the same way, Eqs.(4) and (9), along with 



R R 



B, 



tJV-l 




+Bl 



R R 



(10) 



allow us to rewrite the RHS of the above equation 
as 



(1) iVi?r"i(r„-An)-iV(iV-l)BUii?r"'- (11) 



One important quantity for a condensate made 
of N composite bosons Bq is its normalization fac- 
tor. Let us write it as 



(12) 



If Bq were an elementary boson creation operator, 
we would have Fp^ = 1. For composite bosons, Ff^, 
equal to 1 for TV = 1, decreases when N increases, 
due to what we called "moth eaten effect" [9] : more 
and more free pair states are missing in the op- 
erators of Bq^\Q) due to the Pauli exclusion prin- 
ciple between these A'' pairs as if N little moths 
had eaten these free states. 
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To calculate Fn, we first note that {0\B^bI^\0} 
also reads {0\B^-'^BoBI^\0). We then use Eqs. 
(10, 11). For TO = 1, we find 



F. 



N 



F, 



N-l 



{N-2)\ 



(13) 

and we iterate using Eqs. (10, 11). This shows that 

the Ftv's are linked by 

Fn = i^Ar-i-(iV-l)TiF^_2 

+{N - 1){N - 2)t2Fn-3 + ■ ■ ■ 

+{-lf-\N-l)\rN-iFo (14) 

Eq.(13) also shows that F^ is a decreasing func- 
tion of A^: The moth-eaten effect gets larger and 
larger when N increases. Indeed, the last matrix 
element is positive as seen by expanding it on free 
pair operators; this matrix element then reads 



E 



the sum being taken over different (ki, • 
due to the Pauli exclusion principle. 



(15) 



As a result, the (/)k distribution of the correlated 
pair operator defined in Eq.(16) is related to 

the mean value (N) of the number of up or down 

spin electrons N = J2k '^kf^'^l" I'pBCs) state 

through 



+ 



(21) 



Turning to the fluctuation of this mean value, 
we find that it reads 



(22) 



Since each sum over k is proportional to the sample 
volume, the above ratio goes to zero in the large 
sample limit as 1 over the volume. So, this fluc- 
tiiation is indeed very small which proves that the 
N distribution in the \ipBCs) state is very much 
peaked on its mean value (N). 

We can also consider the mean value of the two- 
pair operator /3k/3k' \4'bcs) state. For k ^ 
k', we find 



III. BCS ANSATZ 



BCS) 



Mk^kUk'Wk' 



(23) 



Let us introduce the unnormalized correlated 
pair creation operator 



E'^p/^^ 
p 



(16) 



with (0|CC^|0) possibly different from 1. From it, 
we construct the following linear combination of 
A''-pair states[18-20] 



+ 00 

E M^^^io) 

N=l 



(17) 



This state also reads 

|<^)=e^>)=np(l + </,p/3t)|0). (18) 

since = due to the Pauli exclusion principle. 
By writing (j>p as Vp/up with \Up\ + \Vp\ = 1, we 
get the usual form of the normalized BCS state as 
product of p operators 

\<l>BCs) = 7l<^) = IlK + «p/5p)|0)' (19) 



where the normalization factor reads 7 = 

Using this p product, the k electron distribution 
in the \^bcs) state is easy to find as 



(Nk) = {4'BCs\at-fau^\4'BCs) 



= \vl\ = 



1 + i^^i 



(20) 



where Fk, often called "pair wave function", is de- 
fined as 



Fk = Uk^k = 



1 + 1 



(24) 



IV. DIRECT CALCULATION OF THE 
A^-PAIR STATE PROBABILITY IN THE 
BCS ANSATZ 

The above calculation of the {N) fluctuation us- 
ing the grand canonical form \(pBCs) of the BCS 
ansatz as a p product is quite convincing to con- 
clude that the A'^-pair state distribution in this 
ansatz is very miich peaked on the N mean value. 
However, a direct calculation of the probability dis- 
tribution, using the N sum is of interest to 
understand the microscopic origin of this peaked 
value. 

If we only consider the {1/N\) prefactor in the 
1^) sum, we could naively conclude that the = 1 
state dominates the sum. However, in order for the 
prefactors of the A''-pair state in the \(p) expansion 
to have some physical meaning, this N-pah state 
must be normalized. 

To do it, we first introduce the normalized cor- 
related pair operator 



(25) 
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where ipp = (jyp/VS and S = \(t^p\, in order to 

have {0\BoBI\0) = 1. The normalized A''-pair state 
associated to the correlated pair operator then 
reads 



iV'iv) 



^0 



|0) 



(26) 



with Fn defined as in Eq.(12), in order to have 

The above equation allows us to rewrite \(pBCs) 
given in Eq.(19) as 



\(pBCs) 



^2. 

N 



FnS^ 

m 



N 



We then note that 



(27) 
since 



{4'BCs\<f>BCs) = 1 and {iJjNli^N) = 1; so, |A^| in- 
deed is the probability distribution of the TV-pair 
states in the BCS ansatz. 

To show that this N distribution is peaked, we 
first note that 1/A''! and Fn both decrease from 
1 when TV increases. So, the ratio F^/Nl also 
decreases when A'' increases. In order to show 
that this decrease is compensated by the increase 
of 5^ in order for A^v to be possibly peaked, we 



first note that the sum S = 



than J2r 



(N) since |Mp| 



I is larger 
1 so \ul\ 



creases with A^. This 5* increase first dominates 
the decrease since, duo to the Stirling for- 

mula, S^/m - {S/N}^ '. So, in the absence of the 
Fjv factor, the Ajv probability distribution distri- 
bution would be peaked on A''** ~ S, which is far 
larger than the pair number mean value (N). 

The moth-eaten effect induced by Pauli block- 
ing on Cooper pairs is responsible for bringing the 
Xn peak back to (N), through the Fn decrease it 
induces. To show it, we first note that this F/v 
decrease does not affect the A^r behavior so much 
for small since, as seen from Eq.(13), Fm /Fm-i 
stays close to 1 for small A^. By contrast, Fm is go- 
ing to play a key role for largo A^'s by changing the 
peak of the probability distribution from A''** = S 
to A''* which, as a definition of the Ajv peak, must 
be such that Ajv»-i — Xn*- Eq.(27) then gives 



Fn* S 



bN* 



(28) 



If Cooper pairs were elementary bosons, Fn would 
stay equal to 1 for all A^ and the peak of the Xn dis- 
tribution would take place for a pair number equal 
to S which is far larger than (JV). For composite 
bosons, the Fn/Fn-i ratio is smaller than 1 since 
Fn decreases with A^, as previously shown; so, Ajv 
is peaked on a A^ value smaller than S. 

If the N* peak were in the dilute regime, the 
Fjv ratio would be close to 1 and A''* would be 



close to S ^ {N). This shows that the solution 
of Eq.(28) is in the dense regime, with Fn/Fn-\ 
substantially smaller than 1. This has to be con- 
trasted to excitons for which the Fn/Fn-i ratios 
always are close to 1: at large density, excitons 
would dissociate into an electron- hole plasma[2]. 

More demanding is to go further and to precisely 
relate the pair number A^* corresponding to the 
Ajv maximum to the pair number mean value {N) 
calculated within the p product form of the BCS 
ansatz. With this goal in mind, let us first cal- 
culate the k electron state population in the B^^ 
condensate. 



V. k-ELECTRON POPULATION IN 
iV-PAIR STATE 

The k electron distribution is straightforward to 

obtain in the \(j)BCs) state; the result is given in 
Eq.(20). To calculate the k electron distribution 
in a A^-pair state iV'jv) is not as easy. 

As a first idea, we could perform a brute force 
calculation of 



(29) 



in the normalized state \iPn) given in Eq.(26), us- 
ing the commutator formalism developed in section 
II. To do it, we start with [ok, Bq\ = v'k&^^k which 
gives by iteration 



[ak,sH = iV<Pk6Lk5r"' • 
This allows us to rewrite Eq.(29) as 



(30) 



N\F, 



N 



(31) 



1 



,,6_k and iterate the 



We then use b-\J}^_-^^ — x — 
process. This gives {N\^n through the following 
Fjv expansion 



N 



\^l\FN-l-{N-l)\^i\FN-2 



N 
Fn 

+{N-l){N-2)\^i\FN-3- 



.(32) 



Eq.(14) shows that the sum ovct k of the above 
bracket reduces to Fjv; so we do have 



$^(^k)jv = N, 



(33) 



as expected. However, through this Fjv expansion, 
it is not easy to show that (Nu) n indeed reduces 
to (Ak) for A^ equal to the peak value A'*. 

A better way to make such a link is to follow 
Leggett [20] and to introduce the operator de- 
fined in Eq.(16) with the k state removed from the 
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sum, namely. 



cl 



(34) 



We then construct its normalized form B]^ = 
Ck/\/5k with = X^p^k l?^pl order to have 
(0|i3k-B^|0) = 1. The associated A^-pair normal- 
ized state then reads 



IV'iV.k) 



4^10) 



ct^io) 



(35) 



where, as for Ffq, the normalization factor .Fjv.k 
is defined as {Q\B^B\^\0) = N\Fn,u in order to 
have {iPnm\'>Pnm) = 1- 

By writing C^^ as (C^ + (/'k/Sk)^ ^^'^ noting 
that (/3k)^ = due to the Pauli exclusion princi- 
ple, wc easily find that the normalized A''-pair state 
\tpN) defined in Eq.(26) also roads 



IV'iv) = 



\4>N,u) + XN,k<l>ul3l\-4'N-l,k) 



N,M 



(36) 



with a;jv,k defined as xn va. Eq. (28), namely 



^iV,k 



fjv-i,k N_ 



(37) 



Using this new expression of {iPn), it becomes 
easy to write the population of the k electron state 
in the A'' pair condensate in a compact form as 



(A^k)iv = 



1 



"W,kl 



(38) 



When compared to (iVk) calculated in the \(j)BCs) 
state, as given in Eq.(20), we see that (-/Vk)Ar cal- 
culated in a A'-pair condensate reduces to (iVk) 
for a;|fk = 1- We then note that, for large sam- 
ples, the number of k states is very large; so, 
the S sum does not change very much when one 
state is removed, i.e., S ~ S^. In the same way 
Fjv — -Fiv.k- As a result, the condition a;^k = 1 
also reads 1 ~ = NFn^i/SFn- This is fulfilled 
for TV = TV* defined in Eq.(28). This shows that 
(Ak),v = (A'k): as expected, the k-electron state 
population calculated in a TV-pair state is equal 
to the one calculated in the grand canonical state 
\<t>BCs) provided that TV corresponds to the max- 
imum value TV* of the TV-pair distribution in the 
BCS St ate. 

To get a better understanding of the link which 
exists between \(j)BCs) and its TV-pair state projec- 
tion I* at), we can also calculate the mean value of 
the two-pair operator /^k/^k' in this |^'jv) state. As 
for ak^akf, a brute force calculation of this mean 



value would give it as a Fjv expansion, not easy 
to compare with Eqs.(23, 24). We can instead fol- 
low Lcggctt [20] and introduce the operator in 
which the k and k' states are missing, namely. 



P5^(k,k') 



(39) 



We again construct its normalized form -Bkk' ^ 

Ckk'/V'S'kk' where 5kk' = Z)p^(k,k') l</'pl and the 
associated TV-pair normalized state 



IV'JV.kk') 



^/N\FN,k, 



c 



kk'l 



iV!i^w,kk'^^k' 



(40) 



with Fw,kk' such that (0|B^k'-Slk'|0) = A^!i^Ar,kk'. 

From C'^ = C]^, + (pu^l^ + i^k' /?k' > i* i^ then easy 
to show that the normalized TV-pair state \tpN) de- 
fined in Eq.(26) also reads 



1 



|V'iV,kk') 

+XNMi^'{4>kPl + </'k',5k')IV'Af-l,kk') 

+a;Ar,kk' a:; jv- 1 ,kk' <t^k4>]<i' /3l^l> \tpN-2,kw) (41) 



where X7v,kk' has a similar form as .Tjv k with now 
two states excluded instead of one, namely. 



F, 



Af-l,kk' 



TV 



^AT.kk' — 



F 



N.kk' 



kk' 



(42) 



In order to still have (f/'jvlV'iv) = 1) the normaliza- 
tion factor J\f in Eq.(41) must be equal to 



+<kk'a;:^-i,kk'l<^kll'/'k'l -(43) 



1/2 



If we now note that, for TV large, a;jv,kk' and 
XN-ikk' are very close, this normalization fac- 

1/2 



tor reduces to a product 



XN,kk'\(l>k'\ 



1/2 



1 + XN,kk'\(l>i\ 



1 + 



It is then easy to show, using 



Eq.(41) for |V'jv), that 
(^iv|/3k4,|Vjv) = 



XNMk"Pk 



XN,kk' 0k' 



^ + X%,kk'\^l\ 1+3^' 



JV,kk'l'/'k'l 



(44) 



When compared to the mean value of /3k/^k' 
culated in the |0bcs) state, as given in Eqs.(23, 
24) , we find that these two results are identical pro- 



vided that 1 = 



Af.kk' 



We then note that x 



N,kk' 



and xjf again are very close because the missing 
(k, k') states do not very much affect F^ and the 
S sum in the large sample limite when the number 



6 



of k's is very large; so F]v,kk' — Fn and ^kk' — S. 
We thus end with the same Pj^PI^, mean values in 
the N pair state IV'at) and in the {(pBCs) state pro- 
vided that A'' is equal to the maximum A'^* of the 
Ajv distribution of A''-pair states in the BCS ansatz, 
which also is the mean value of the particle number 
(TV) in the \4>bcs) state. 



VI. COOPER PAIR WAVE FUNCTION 

We wish to end this work by reconsidering what 
should be called "pair- wave function" , since Fk de- 
fined in Eq. (24) is often called this way. Cooper 
pairs like cxcitons arc composite bosons made of 
two fcrmions. Their creation operators thus read 
as a sum of free fermion pair creation operators 
/3^. To properly identify what should be called 
"Cooper pair wave function" , it is enlightening to 
compare their creation operator with the exciton 
creation operator. 

We commonly distinguish two types of excitons: 
Wannicr cxcitons[21] and Frcnkcl cxcitons[22, 23], 
the latter having closer similarity with Cooper 
pairs. Indeed, Wannier excitons are constructed 
on free electrons and free holes, so that they arc 
"double index" composite bosons. Their creation 
operators read 



(45) 



with i = {Qi, Vi), where Qi is the exciton center-of- 
mass momentum and v-i its relative motion index. 
The prefactor (k/j, ke|«) of the i exciton expansion 
on free electron-hole pairs, is the i exciton wave 
function in momentum space. 

By contrast, Frenkel excitons arc made of atomic 
excitations for atoms on a regular lattice. They 
thus are "single index" composite bosons, their cre- 
ation operators reading as 



(46) 



the 



R„ is the position of the excited atom and Ng 
number of atomic sites. 

Cooper pairs also arc "single index" composite 
bosons since an up spin electron k is coupled to one 
down spin electron (— k) only by the reduced BCS 
potential. The (normalized) creation operator of 
one Cooper pair in the BCS condensate reads, as 
quoted above. 



B, 



-k' 



(47) 



where ip\^ = <?i>k/v^ with 



^'k = and 5 



It view of Eqs. (45-47), </5k must be taken as the 
Cooper pair wave function, without any ambiguity, 
in spite of other quantities like _Fk ~ uJ^Uk often 
quoted as "pair wave function" in the literature. 

Although Cooper pairs and Frenkel excitons 
both are "single index" composite bosons, they 
however have a few major differences. One of them 
comes from the fact that the Frenkel exciton wave 
function is just a phase, so that the distribution of 
the excited site n in a Frenkel exciton is flat. By 
contrast, the Cooper pair wave function is given 
by ^pk = (l)k/ VS with 



3k = 



Vk 
Uk 



A 



Ek + Ck Ek + ^k 



(48) 



= — ^ and = + A^, the electron k 
energy being Ck = k^/2TO while /z is the chemical 
potential of the \(f>BCs) state in the grand canoni- 
cal ensemble. In the usual BCS configuration with 
a reduced BCS potential extending symmetrically 
on a phonon energy scale on both sides of the nor- 
mal electron Fermi sea, /x is in the middle of the 
layer over which the potential acts, this layer ex- 
tension ft being twice the phonon energy (see Fig. 
2). The gap A then reads in the small potential 
limit as A ~ Qe~^^''°^ , where po is the density of 
states taken as constant over the layer where the 
potential acts. As a result, ^k has three different 
scales: 0k — 1 for Ck very close to fi on the A scale. 
(j)k — 2e^/P°^ {pL — ek)/fi for ek far below /i o n the 
A scale and l/(/>k ^ 2e^P°^ {e-^ - for ek far 

above /x on this scale. 




FIG. 2: The reduced BCS potential extends over an 
energy Q, of the order of twice the phonon energy, above 
a non-interacting electron Fermi sea with Fermi energy 
cfq. The dashed line indicates the chemical potential 



If we now turn to the normalized distribution, 
i.e., the Cooper pair wave function in momentum 
space, it reads v'k = (pk/VS where, using Eq. (48) 



^ = Ei^i 

k 



iVn(l- 



6A2 



) 



,2/poV 



(49) 
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= Po^ is the total number of pair states in the 
potential layer for a constant density of states po- 
Since e~^^P"^ is very small in the small coupling 
limit, the Cooper pair wave function is sizeable 
between ep^ and — A/2 only, where ep^ is the 
Fermi energy of the non-interacting electrons (sec 
Fig. 2). (^k then scales, within irrelevant numerical 
prefactors, as 



1 M - Ek 



(50) 



It has a small tail (see Fig. 3) of the order of 



(2) e" 



(51) 



for electron energies in the range ±A around /z. 

For higher energy, i.e., for ek between /i + A and 
epo + ^, the wave function is even smaller, being 
of the order of 



(3) 



VNq ek - m ' 



(52) 



This shows that the sizeable part of the Cooper 
pair wave function 93k, which is the normalized 
form of 0k = Wk/wk, is a linearly decreasing func- 
tion of ek between the non-interacting electron 
Fermi sea ep^ and the normal electron Fermi sea 
ep = €p„ + n/2, in the case of a potential extend- 
ing symmetrically on both sides of this Fermi sea. 
The number of pair states with sizeable weight in 
the Cooper pair distribution thus is of the order of 
Nq/2. This understanding has to be contrasted 
with what is often called " pair wave function" , 
namely _Fk = ^^k^^ki £^nd which is highly peaked at 
€p (see insert of Fig. 3). i^k is physically related 
to the excitation of electron-hole pairs in the BCS 
condensate while tpk is associated to the ground 
state of up and down spin electrons added to the 
non- interacting Fermi sea ep^ and paired by the 
reduced BCS potential. 

For completeness, in addition to Cooper pairs 
with creation operator making the BCS con- 
densate, we must mention the "single Cooper pair" 
creation operator derived by Cooper when study- 
ing a single pair of up and down spin electrons 
added to the ep„ Fermi sea. Its ( unnormalized ) 
creation operator reads 



E 



2ek - El 



-ay 



k- 



(53) 



the sum being taken over ep^ < ek < e^o + The 
single pair binding energy is Ei = 2epg — Cc, where 
ec 20e-2/^o^ for smaU V. The above equation 



shows that the wave function of the i?]v=i|0) state 
is concentrated on a ec scale above ep„\ so, that 
the amount of pair states with sizeable weight in 
-Sjv=i is of the order of = poCc which is far 
smaller than the number of pairs in the Bq opera- 
tor making the BCS condensate. 




FIG. 3: The Cooper pair wave function ipk essentially 
has a linear decrease between sfq and efq + § — A* and 
a very small tail between fj, and epg +Q. 
Insert: the "pair wave function" defined as -Fk = Mk^k 
is concentrated on a A scale around sfq + §■ 



VII. CONCLUSION 



In this work, we address the BCS condensate of 
Cooper pairs in the canonical ensemble. To eas- 
ily handle the Pauli exclusion principle between a 
given number of paired electrons, we first develop 
a formalism appropriate to Cooper pairs which is 
based on a set of commutators. We then use it 
to, in particular, show that the Pauli exclusion 
principle between Cooper pairs is fully responsi- 
ble for the correct value of the probability distri- 
bution peak for iV-pair states in the BCS wave 
function ansatz. The standard grand canonical en- 
semble approach tends to mask the key role played 
by Pauli blocking in this problem. We end by re- 
considering what should be called "pair wave func- 
tion" through a comparison with other composite 
bosons like Wannier and Frenkel excitons. 
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